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I. INTRODUCTION 


-o- 


Combat active replacement factors, on GAkrs are Loewe 


tics planning factors currently used by the U5 Marine Gore. 


as estimates of equipment losses in future conflicts. Their 
values have a Significant impact on procurement, stock- 
piling, and plans for shipping requirements. A replacement 


ut 


factor is defined as the estimated percentage of equipment 
in use that will require replacement during a given period 
due to wear-out beyond repair, enemy action, abandonment, 
pilferage, and other causes, except catastrophies”" [Ref. 1]. 
Further, combat active replacement factors are applied for 
units during those periods when they are actually in active 
combat operations. A force in Gontact with the « emenc samc 
considered to be active combat. 


There are alternatives to the use of computerized war 


games and simulation to obtain values for CARFs. [It is 
possible to produce estimates directly, perhaps using 
professional military judgement and experience. In some 


cases CARF values so generated may be preferred to those 
obtained from combat models, since there may be more clarity 
about what the number was based upon, and what considera- 
tions went into its estimation [Ref. 2]. 


A variation on estimating CARFs directly is to use 


professional judgement to estimate a related measure, the 
mean time until loss or MTTL. The MTTL is the average time 
one would expect an item to survive in the combat 
environment. 


This thesis discusses some ways of generating CARF 
values without to resort to war games or to simulation or 
direct estimation. Basic to this method 1S the estimacionm 


of values for MTTL, the mean time until loss, or the 


estimation of the intensity function for a nonhomogeneous 
Poisson process, andthe method will be most applicable to 
Situations where the MITL can be estimated with more confi- 
dence than the related CARF can be directly estimated. 


There are three major types of life distributions that a 


certain item might follow in the real world. These distri- 
butions are the exponential, the Weibull, and the gamma. 
First, we will look at the loss process characteristics of 


life distributions and then derive the CARF expressions for 
several kinds of scenarios according to each life distribu- 
tion. These Scenarios by which losses could occur in combat 
are: 


Ll. Cases where all items of that type are vulnerable 


icidawieeand throughout the wveombat period, at the 
same MTTL values, 

2. Cases where subsets of the items in use _ have 
different MTTL values, 

3. Cases where the MTITL changes at some designated time 


during the combat period, and 
4. Cases where one item is vulnerable initially and 
throughout the combat period at the same MTTL, but 
replacements are not vulnerable until put into uSe. 
Regardless of the life distribution, we will also 
discuss the case where an item's loss process follows a 
nonhomogeneous Poisson process at a loss rate that is time 
dependent. 
imethem.Lollowing chapter, “we will discuss™loss protCess 
Beoreeteristives for each Wetewdastributiony andwftor a nonho- 


mogeneous Poisson process. 


ii; LOSS PROCESS sCHARAGTERISi Gs 
Differences in the loss process may be portrayed by 
several kinds of life distributions. The loss rate is 
constant for the exponential, but for the Weibull and gamma 
distributions, changing loss rates can be represented by the 
distribution parameters. In this chapter we will look first 
at some loss process characteristics for each life distribu- 


tion and then for a nonhomogeneous Poisson process. 


A. LOSS PROCESS FOR THE EXPONENTIAL 


The exponential distribution is in several senses the 
most fundamental distribution in reliability theory. We 
will look at a loss process where the event that an item is 
lost is independent of the loss of any other on-line item of 
the same kind, and where the chance of a surviving item 
being lost on the next day is independent of how many days 
it has already been in combat. This means that the indi- 
vidual item's loss rate may be considered constant over the 
period of time we are examining. For constant loss rates 
and independent losses, the probability distribution for the 
time an item survives in combat will be the exponential 
distribution (Rete 24. The expected value of this random 
variable is the item's MTTL. Let T be the combat survival 
time for a specific item, and let R be the scale parameter 
for the exponential. (For the exponential distribution, R 
is also the loss rate, with units of items per day.) The 
density function for the time t until the item is lost is 
exponential: 


=INE 
F(t) = Re R20, t>0 


LO 


From this we may immediately obtain the probability that the 
menets lost on or before time t, F(t), as 


= Nee 
EGpm= p( lsmp- =e ; 20) 


~ - 


We also have 1 - F(t), the probability that the item is 
still surviving just after time t, as 
ie Ge ep(T>t) = ia 
Since the expected value or average value of an exponen- 
tially distributed random variable is the reciprocal of its 


parameter, we now have the mean time until loss MTTL as 
Me biee = lh ) KR 


meecordingly, the probability that the item is not lost 


aourane D days of combat is 


=D) MrerL (2 a) 
=e : 
In the following section, we W1lll look at some loss 
process characteristics for thie Weibull ‘ and gamma 
Seestributions. 


Pree cess PROCESS FOR OTHER LIFE DISTRIBUTIONS 


Other parametric families of life distributions arising 
in combat situations may be constructed by assuming a loss 
rate as decreasing or increasing rather than being constant. 

In each case, the loss rate 1s considered monotone over 
the time. For decreasing or increasing loss rates, useful 
meobability dastributions for the time an item survives in 
combat are the Weibull or gamma . We will look first at the 
weibuliadistribution. 


JE 


l. Loss Process fonvtnes wer pu 


The loss rate r(t) for the Weibull distribution wien 


shape parameter C.and scale parameter R is 


A-J - 

r(t) =seee se) MW >0, R20, t>0 , 
and it iS imecreasime "for ce decreasing for 0<C<l)) lame 
Constant for Comper ese The density function for the 


time t€ tntid “the wrtem 1S host 25 


vo. 
A A-1 -(Rt) 
le e 


£() oe t20 


From this we may obtain the probability that the item is 


lost on or before time E5u Cees 


eG 
SME 
F(t jc =3p Crt p= sae : ee) 
We also have l - F(t), the probability that the item is 


still survivine@itist arter time taweas 


1 - F(t) = as 


We have the mean time until loss MTTL as [Ref. 4] 


Accordingly, if time until loss is distributed according gee 
the Weibull distribution, the probability that the item is 


not lost dusine > days of combacet- 


o& MTITL (2.5 


ieZ 


Maes is equivalent to (2.1) for the exponential when & is 
ieesance [°(1) = 1. : 


Vememoss Process for the Gamma 


|e Se 


oad 


The loss rate r(t) for the gamma distribution with 
shape parameter A and scale parameter R is increasing 
meme > 1, decreasing for 0 <A< 1, and constant for &= 1 


[Ref. 3]. The density function for the time t until the 
item is lost is 


f(t) = ------------ : Se OeR= 0, 120 


Meee probability that the item is lost on or before time t, 
Eat), ks 


ONCE eS os a ae, t20 


When “A is a positive integer, F(t) may be written in closed 
merm as [Ref. 3] 


ni 
a JRE = |e 
AR es a ey e ; ee8, 


Mieessurvival probability that the item is still surviving 


Mmlist after t, is 


and from the above formula, the scale parameter, R, can be 


derived as 


i 


According to Ehisy the probaieeirc yee ie ee ieeecn ceo t vce 
during D days of combat is 


sl 
(A.D/MTTL) -(a& /MTTL)D (2 25 
DP. = 22 eee es e 

T= i! 
When Coase lee nis is equiva lenemmeton (2.15) for the 


exponential. 


C. NONHOMOGENEOUS POISSON PROCESS OF THE LOSS 


In this section we consider the nonhomogeneous, also 
called the nonstationary, Poisson process, which is obtained 
by allowing the loss rate to be a function of t. As a 


prelude to giving a definition of a nonhomogeneous Poisson 


process, we shall define the concept of a function f(.) 
being O(h). The funetion £(.) is Said to be O(h) [Re£t 2m 
alate 
fla 
lis, Speen 
h->0 h 
inorder tor’ ene fumctton. &( 2) to be O(h) it is necessary 


that f(h)/h go to zero as h goes to zero. But if h goes to 
zero, the only way for f(h)/h t6 go to zero is for £({h )teeme 
go to zero faster than h does. That is, for h small, "fG 
must be small compared to h. 

According to the definition of Ross [Ref. 5], the 
counting process {N(t),t20} is said to be a nonhomogeneous 


PoussonM process: Wert Int enor yet Um ero nee Geo) Ont S 


Le: “NCO, =20R 
2. {N(t),t20} has independent increments. 
35° pt hy G2 j= On 
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pe (Ceti j-N(t)sme} = X(t)h + O(h). 
If we let | 


me) = [iacsyas | 


Maem tt can be shown that 


NeCeE Pen” Sam 
ee ee et? Cees Pmt) ed les. 


, 
In other words, N(tt+s)-N(t) is Poisson distributed with mean 
Meets -m(t). Tats, N(t) issPoisson distributed with m(t), 
and for this reason m(t) is called the mean value function 
of the process. Note that if x = 0, there are no losses 
Prom t to t+s, and thus p{N(t+s)-N(t)=0} is the probability 


that one item survives for this period, and we have 
p(item survives from t to t+s, given survival to t) 


-(m(tt+s)-m(t)) (2 a5) 
=o : 
With these characterizations of the above general loss 
processes established, we will ocwrat “wanmvous locs 
Scenarios to derive CARFs for the exponential in the 


following chapter. 


if 


IIIT. CARF GENERATION MODEL FOR THE EXPONENTIAL LIFE 


Sr ro LS se 


DISTRIBUTION oe 


The case where the life distribution for an item is 
exponential has been examined by Lindsay [Ref. 2), Wit 
combat active replacement factor expressions developed for 
several scenarios. 


In this chapter we will summarize his results. 


A. CARFS WHEN TTEMS HAVE THE SAMESIa TU 


We look first at the case where all items are committed 
and initially vulnerable at the same MTTL throughout’ the 
D-day period. If n items are the initial in-use amount and 
are committed with independent losses, then the number of 
items x that would be lost during the D-day period will be 
binomially distributed, and the average number (out of n) 
lost in D days is simply the mean of the binomial dicstimieas 
EHOn. One 

-D/MitE 


Average number lost = np = n(l - e 


From this we may readily obtain a CARF value: 


(in-use amount) 
one 
-D/MTTL (3785 
GARB = “(jet ae er OeD 


Another use of (3.1) would ask if the existing CARF 


value would yield an MTTL that seemed reasonable, using 


MTTL = -D/1n(l - CARF/100) 
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The association between CARF and MITL for this case is 
illustrated by the values in the exponential column of 
Appendix A. 


-— 


Peek ro WHERE NOT ALL ITEMS HAVE THE SAME MTTL 


If all items are Paola lly oviilnerable, bute awasen 
different MTTL values for the D-day period, then CARF gener- 
ation from the MTTL values is an extension of the case where 
all items have the same MTTL. 

fmeeeOpoEEions. Pl, pws, pk Of the in-use amount; n, 
mewesmean times until loss MITL1, MITL2, ... , MTTLk. The 
average number lost in D days for any subgroup 1 1S 


= by ME lioY 
i lel - ©& ) 3 
Al 


ema thus 


k 
: y -D/MTTLi (352) 
CAR Cee (1 = 4 Cue) —~: 
i=l 1 

This is, of course, simply a weighted average of CARF values 
from the previous caSe. An example of the CARF values for 


this case is shown in the exponential column of Appendix C. 


weeeeak?S WHERE THERE IS A CHANGE IN THE MTTL 


It may Sometimes be of interest to construct a CARF for 
a Situation where there is a change in the combat scenario 
Gunmeme the D-day period. One example of this is the case 
where the first portion of the D days is an amphibious oper- 
ation, and the MTTL might subsequently change. 

fem this type of situation, we let 1/MTTL1"be the loss 


rate for the first Dl days where xl items are lost, and 
Sei? Sbewsthe loss rates foruthe remaining (D-Dl) days 
where x2 items are lost. Here, total losses for the D-day 


merited are xl+x2. 
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Let ql be the probability that am item woula=™ be” lose 
during the first period Dl and q2 be the probability thar am 
item would be lost during the remaining D-Dl days. Then 


oo 


Se Ml ene 
ql=l-e 
and 
=( DED 1) / ME riEZ 
qZ2 =l-e 
Both xl and x2 are binomially distributed, xl with parame- 
ters. (nee), and x2 with parameters for the conditional 


binomial distribution) ((n-x! \fa2) eanen 


n 
EL Xi +x2Z oe ale do (nex) eye jelGeil) 
XLl= 


mo q+ GZ (neces) 


(33°) 
=enil = (1 -qi)( hse. 
and the CARF would be 
(Di /MTTL1 = (pep te (3.4) 


CARF = (l-e ye Clero: 


The association between CARF and MTTL for this case is 
illustrated by the “values in the exponential column @eE 


Appendix D. 


D. CARFS WHERE AN ON-LINE ITEM 15 REPLACED BY A PREVIOUS 
INVULNERABLE ITEM 


The expressions for generating CARFS in the previous 
Sections have all been based upon a Situation where all 
items were initially vulnerable. A different Ioss process 
would place one item on line, and structure its (possibly 
repeated) replacement with heretofore invulnerable items. 

If the supply of reserve items is very large, or large 


enough so that the chance of its exhaustion is negligible, 


iS 


Pemedm Structure this case by noting that the loss process 
is Simply a sequence of exponentially distributed time 
intervals over a D-day period. If we ignore boundary condi- 
Eons , the number of such intervals would be a Poisson 
G@eestributed random variable with a mean of D/MTITL, and the 

CARF [Ref. 2] would be 
Davis Cae, C3250) 

CARF = ( ------------- Geos) 
> i Sie eaoriene 

Equation (3.5) may be a reasonable approximation for an 
item with a low CARF value. If the chance of running out of 
replacements is not negligible, however, then (3.5) will 
yield an overstated CARF value. This can be converted (at 
the cost of simplicity) as follows. Let n be the initial 
amount of an item, and x be the losses in D days. Then, with 


the Poisson probability distribution we can write 


x =D/MEEE 
(D/MTTL) e 
a eee ; = tee Oras Agta 
a. 
me(x) = 
x ) (B26) 
7) x =D/ MILE 
y (D/MTTL) e 
aS 9 = Se a HS HS Sees oe = 1) : 
x=n a! 
Mometche loss distribution. Taking expected values, Eas 
yields a CARF expression 
n- oD Mt LL 
100 7 by Meee 
CARF = (---)\| QQ. ------------------- 
n x=0 ee 
gtk (D/MTTL)* -D/MTTL (3.7) 
e 
el (Ib = D re ) : 
x=Q , 


Examples of CARF values for each case are shown in Table l. 
Curves associated with Table 1 are displayed in Figure 3.1. 
It can be seen that the approximate and correct CARF values 


are the same when the initial amount of an item, N, is 10. 


ie 


iAb Lee 
EXAMPLES SHOWING THE IMPACT ON CARF VALUES OF ASSUMING 
INFINITE SPARES FOR A HOMOGENEOUS POISSON PROCESS 
CARE wvaiues, ) 
Mini (davis) correct 
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In the next chapter, we will look at the CARF generation 


model for the Weibull and gamma life distributions. 
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(D = 30 DAYS, N: NO OF ITEMS) 


LEGEND 
APPROX 
COiKmec!) ~~ 


15 
MTTLS (DAYS) 





Eo BU Spee e The Impact _of Assuming 
nfinite Spares on CARF Values. 


Zu 


IV. CARF GENERATION MODED POR O2THER ELVES teens 


Pap ey Spm me me mmm 


~~ 


In this chapter, we will derive CARF expressions for 
several scenarios using the Weibull distribution and the 


gamma distribution. 


A. MODEL FOR THE WELSULL) DISTRIisU rah 


We are interested in looking at the Weibull distribution 
because it is one of the most commonly used life distribu- 


Cilons in rebiabpwbie,  eneon- 


lie CARFs When Items Have the Same MTTL 


When all items are committed and initially vulner- 
able with the same MTTL throughout the D-day period, the 
probability that an item is lost during the D days 1s, “fueq 


Or a 
a 
: ) 
awake 
l-e 2 
Therefore, 


A 


are 
ee (4.1) 
GARE = lee? e 


) (100 ) 


This is equivalent to (3.1) for the exponential life distri- 
DUELoOn. whence Gis wen Equatilon (+e) permits a CARF tomes 
computed from an MTTL estimate for the case where each item 
in the in-use amount 1S committed initially and at the same 
MTTL, and where losses are independent. Equation (4.1) cam 
also be used to determine if the existing CARF value would 


yleld an MTTL that seemed reasonable. Here we could use 


Ze 


id Oe eee D [i (1/A) re aaa 


o& [-1n(l- CARF/100)]%A 
The association between CARF and MTTL for this case 
is illustrated numerically by the values in the Weibull 
columns of the table in Appendix A and Appendix B when co is 
2 and 3, respectively. 


2. CARFs Where Not All Items Have the Same MTTL 


If all items are initially vulnerable, but with 
different MTTL values for the D-day period, then CARF gener- 
ation from the MTITL values is an extension of the case where 
all items have the same MTTL. 


POmmamceneralelLormuvaye Leteproportions pl,pZ,...,pk 
of the in-uSe amount, mn, have mean times until loss MTTLI, 
fees 2... 286,06 | 6MEELK. The average number lost in D.days for 


any subgroup i is 


coy 
a an A MTTLi 
ae 


and we have 


; OetTLi (4.2) 
Con at a 2 - e | yO) CL Oe a 


As an example of the use of (4.2), suppose 20 % of 


Eee Ln-use amount have an MITL of 100 days, half have an 
Mitinot SO days, and the remaining 30 %4 have an MITL of 120 


days. Then from (4.2), 
eerste | oe) 
Z LOO 2-30 
Nites du =e ) 


(22) eee 
; 2 +120 j 


Ger = a0. 241 - 


Oo (aL o= re Jere OO!) s=aro 203) 7, 


when Ais 2 and the combat period D is 30 days. This is, of 
course, Simply a weighted average of CARF “values from tne 
previous case. An example of the CARF values for this case 
is illustrated by the values in the Weibull column of 


Appendix C. 


3. CARFs Where There is a Change in the MTTL 


We next look at the scenario where the Weibull 
distribution for time until loss changes at time Dl, ~~ whaler 
we shall represent as a change in the mean for the Weibull 
distribution from MTITLI to MITLZ. The “structure we 97a 
follow will assume that the loss rate after the mean changes 
to MITL2 begins “with argument DE rather thane: The 
Survival function of a life length T for the Weibull distri- 
bution with shape parameter & and scale parameter R is 

: -(Rt)~ 
F(t) = pC) eee : oh >0, R>O, t20 

Let F;(t) be the probability that an item with loss 

rate E(t) Sumpvivece to tine Ceforeis me... From (2.2)7 eae 


SUrVivVal probabil, F, (D1) that an item survives to time Dl 


1s 


pie a 
; ( ; rete 
e€ ° 
Sanmce R= Coe om ria. The conditional survival prob- 
ability for the time period Dl to D, given survival to oe 
1s 
F,(D) / F (D1) 
A X 
Dice) DL [* (1/2) 
; Bee 7 nate) 
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mhen the probability F(D) BittiMedn meenuesumvyives to D is 


Ee!) = ee F (0) / #001) 


eeMmrril x MITe2 


—— Soe See Ss 


A MITL2 


ey ot. =e 
o METLI A MTTL2 


A 
1 coe C ees D1 ibe) 


=e 


Therefore, 


(4.3) 


T Pe cs i (1 
2 ee aaa Di ee 
AMTTLI AMTTL2 
CARF = \l-e (LOND ie 


When MTTL1 = MITL2, this is equivalent to (4.1) for the 
scenario when items have the same MTTL, and when A= 1, this 
reduces to the result (3.4) for exponentially distributed 
Bombat lives. The association between CARF and MTTL for 


this case is shown in the Weibull column of Appendix D. 


Pee MODEL FOR THE GAMMA DISTRIBUTION 


The next case we shall consider is that when the life of 
an item in combat follows the gamma distribution. We will 


imgg< at three combat Scenarios. 


1. CARFs When Items Have the Same MTTL 


Os SS ae 


This is the case where all items are committed and 
initially vulnerable at the same MTTL throughout the D-day 
Meriod. From (2.3), Ete sprGbabilnme that an. item as: lost 
during the D days is 


ee] i 
y (CSDyiMETL ) -(A/MTTL)D 
=o e 


where a 1S a positive integer. 


a) 


Therefore, 


CARF = (1- 2 (hD/wTTL)” Eel 0 ae 


As with the Weibull distribution, when A= 1 this is equiv- 
alent to (3.1) for the exponential distribution. 

The association between CARF and MTTL for this case 
is illustrated numerically by the values in the gamma 


columns of Appendix A and Appendix B when Ais 2 and 3, 
respectively. 


2. CARFs Where Not All Items Have the Same MTTL 


This is the case where not all items have the same 
MTITL » for “the @)-day = combat spend. If all items are 
initially vulnerable, but with different MTTL values for the 
D-day period, then CARF generation from the MITL values is 


an extension of the case where all items have the same Miia 


For a general formula, let proportions pl,pZ,... 3a 
of the in-use amount, n, have mean times until loss MTTLI, 
Mi TiZ eee , MTTLk. The average number lost in D days for 


any subgroup 1 1S 


eons, (AD/MITLL)” “(A/MITLE)D 
= j! 
and we have 
k A -] (A D/MTTLi \ -(A/MTTLi)D (4.53 
cape 2K » PeGle:= 2. appeals pgs i ve 
= j= 5! 


This is the same as (3.2) for the e€xponentialeuhen ouwsle 
As an example of the use of (4.5), suppose 20 % of 
the in-use amount have an MTTL of 100 days, half have an 


MTTL of 80 days, and the remaining 30 % have an MTTL of 120 
days. Then 


Ae 


PT PICOMGG2 ( Tones, <.-----2----- " 


((2» 30)/100)° (2/100) (30), 
Se af 


2 


op eae 2D ------------ e 
= ae 
: Ge 30)/120)> CZ ZO C30) 
7 10% 8 2 ------------- e ») 
j=0 ae 
Sees la yaran 
when “is 2 and the combat period D is 30 days. MMoulS, ole 


Simply a weighted average of the CARF values from the 
previous case. 
Other examples of CARF values for this case and 


these proportions are shown in the gamma column of Appendix 
Oe 


Pee CaREs Where Ihere 1s a Chanpe in the MITL 


We will look at the scenario where the gamma distri- 
bution changes at time Dl, which we shall represent asa 
change in the mean for the gamma distribution from MTTLI1 to 
MTTL2. As before, we will assume that the loss rate for the 
second distribution begins with argument D1. The survival 
Mimerion of a life length IT for the gamma distribution with 
Shape parameter A and scale parameter R is 


ees UT 


miere —eiS 2 pogitiwe integer. Let F; (t) be the probability 
that an item with loss rate r;(t) Survives to time t. for 


meme. Erom (2.3) ,eethe survival wrebability F, (om hat 


iaelniem Sumvives to Di can be obtained as 


Jia 


i=0 a 


b) 


Since R = A /MIEPE: As with the Weibull distribution, ene 
probability that an item survives for the time period Dl to 
D, given “Sugvaay a patoe » lees 


F,(D)/ #001) 


A-] i 
7 y (A D/MTTLZ ) =(_ Ag Mire? bp 
=" 7 (ae e ) 


—_—_— == = = = es eee lel 


A-1 i 
5 (AD1/MTTL2) -( A /MTTL2)D1 
( e ) 


IBC i! 
and the probability that an item survives to Dis 
F() = SCD) 7 FeO 
(D) LOD #02) / Fv 


Then the CARF could be obtained from 


A -] i 
1 
ees = ( a ~ (A D1L/MTTLL) p ny MITE) Dp 


1=0 i! 
A-1 at 
.. (G20) / Mitr ) =( ot. / MiEIE2 )D 
eA (cee & lre 7° ) 


a=] i 
/ ee ee ee ) (100) 


1=0 a 
When Miijei Saale 7s this is -the Same CARF expression as 
(4.4) for the scenario when items have the same MTTL, and 
when A= 1, this reduces to the result (3.4) for exponen- 
tially distributed combat Lives. The association between 
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Peewee ti toOretnisecase “15 illustrated by the values in 
the gamma column of Appendix D. 
In the next chapter, we will look at the CARF gener- 


ation model for a nonhomogeneous Poisson process of an item 
loss. 


ao 


V. CARF GENERATION MODEL FOR A NONHOMOGENEOUS POISSON 


[SSS mm mS o—E—E—————————=£&&£&££EeeEI!“7“"“—_ ES 


PROG Oo a= 


Regardless of the life distribution, another case of 


interest is where an item's loss process follows a nonhomo- 


geneous (nonstationary) Poisson process. The intensity 
function, A(t), related to this monhomogeneous  Poisosen 
process isthe Joss rate as 4a funcercnmo: et  [pRetwole We 


will look at four combat Scenarios. 


A. CARFS WHEN ALL ITEMS HAVE THE SAME INTENSITY FUNCT IOR 


We will look first at: the combat Scenario where all 
ltems are committed and initially vulnerable at the same 
intensity function of the nonhomogeneous Poisson process 
throughout the D-day combat period. If mn items are the 
initial in-use amount, and each has a probability ¢ thaceem 
item would be lost during a D-day period, then the number of 
items that would be lost will be binomially distributed, and 
the average number lost in D days is the mean of the bino- 
mial distribution: 

From equation (2.5), the probability q that an item 
would be lost before D days is 

qg =l ag ; (5.1) 
Thus, the average number of items being lost during the 


combat period D is the binomial mean 


n( 1 - a), 
Therefore, 
CARF = (l - a), (190) ae 
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inMvemcasesgummacsacaatedswith this process will be shown 
here according to the shapes of the intensity function. We 


feet rarst look at the linear intensity function, 


Beee)=9a + bt | a20% D=OReIO=<E=D , 
where "a" is the intercept and "b" is the slope. When b = 
QO, this scenario could be treated as the same scenario for 


the exponential since the loss rate becomes constant. 


D(30 days) 





Figure 5.1 Linear Imtensity Function. 
For example, as shown in Figure 5.1, if a = 0, b = 1/150 and 
D = 30 days, then 


eee) = (1/150)t , 
and 


30 
m(30) = I; (1/150)t dt = 3 


eeom this, 
= -3 
CARF = (1 - e )(100) = 95.0 &% 


This example might be applicable in a combat scenario where 
continuous reinforcement forces on the enemy side are 


expected. 


Sil 


A quadratic form of the intensity function  ioethe jeace 


where the intensity function is gradually increasing in 


the form Gian quadratiemeas time goes on and then 
decreasing as Shown in Figure 5.2. This mtght be appli- 
cable in @ general) compaes scenario In this case, the 


intensity function, A(t), can be expressed as 
Z 
(CC) =a er) 2) eee a>Oeee-0, Osta rae 
where D = 24/b/a. The mean value function, m(D), could be 


abeaimed sas 


D Z 
m(D) = [eee (ae= Dy2)) + sbjae 


= -a(t - D/2)*+ 





Devoe. Je Quadratic Form of the Intensity Funeerene 
For example, let the quadratic form of the intensity func- 
tion be 

SOC Ns -1/2250(t-15)°+ O.1 
Then 


30 y) 
m(30) f (21/7 2950(t=15) On bya 


Therefore, 


BZ 


2 
wpre= (1 =e yC100) = 86.5 Z 


PeeeeCARFS WHERE NOT ALL ITEMS HAVE THE SAME INTENSITY 
FUNCTION 


If all items are imieitaliy vulmerable, but with 
different intensity functions for the D-day period, then the 
CARF generation for this scenario iS an extension of the 


case where all items have the same intensity function. 


met proportions pl,p2,..:,pk of the in-use amount, n, 
Meaemintensity functions Al(t),A2(t),...,Ak(t) and mi(t) 
be a mean value function for subgroup i. Then the average 


number lost in D days for any subgroup i is 


-mi(D) 
pi n(l - e 


PS , 


k e 
CARF = om pi(l - _ We) COO rn «. vote) 
L= 


SemeecaAkrS WHERE THERE iS A CHANGE IN THE INTENSITY FUNCTION 


The nonhomogeneous Poisson process approach to CARF 
development offers an opportunity to handle an intensity 
function which is piecewise. For example, if the loss rate 
1s linearly first increasing and then decreasing during the 
Deaay Period, the intensity function is 

ay FWD ts : 0 pee > Oe << DL 

A(t) = ’ t ' 

a eee Dll). geesOme tb <0. D1<t=D 


[Mae mean value functionm can be evaluated as 


CD) = p> + bt)d fe se 5)! Dl))d 
m = ; (Cat Evade ae (t- IGE 


2/6) 






0.015 =. 9 7aaaeet-15) 


D(30 days) 


Figure 5.3 Linearly Increasing, then Decreasing 
Intensity Punceron- 


For example, if the intensity function is 


Oreo uiats : OS. \6 dls 
A(t) = 
Ones -0 0CLtea)) - USS ees 0 


lA 


as shown in Figure 5.3, then 


I 


IES £119) 
men) I CO 00M e dt - [ $°0.015-0.001(¢-15) ae = 0.2228 


From this, 


CARF = (1 - 5 ee (1007) ae Ol 

The case where there is a change in the constant loss 
rate, as shown in Figure 5.4, could also be treated 
nonhomogeneous Poisson process. However, this sinpds, 
provides another way of deriving results given earlier when 


time to loss was exponentially distributed. 
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Figure 5.4 The Constant Loss Rate Changed 
During the Combat Period. 


D. CARFS WHERE AN ON-LINE ITEM IS REPLACED BY A PREVIOUSLY 
INVULNERABLE ITEM 


The scenario where an on-line item with an intensity 
mumetion A(t) Pome placcd Meer mricncmellome invulnerable but 
identical items may also be examined when losses are 
described as a nonhomogeneous Poisson process. A loss 
Process would place one item on line, and structure its 
replacement with previously invulnerable items. in other 
mereeas, when an item is brought in at time t, its initial 
loss rate is assumed to be A(0) rather Shan. » Aes). 

Let n be the initial amount of an item andx _ be the 
losses in D days. iijeene Chaneesereruinning outsol replace- 


Memes is not negligible, then with formula (2.4) for the 


Honmomezeneous Poisson probability distribution, we can 
write 
m(D = OB) 
0) a 
rei 
p(x) = Cue 
¥ oe 7 no) 
3 m(D ) - m(D) 
------ a ‘ xi on ; 
X=n xm 
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where 


D 
m(D) = [ oxerat 
In a way similar to that done for the homogeneous 
Poisson process discussed in Chapter III, a CARF expression 


can be obtained by taking the expected values as 


100 a (D)” (D) 
x Mm -m 
CARF = (---) eee e 
n x=0 me (5 769) 
Se m(D)* = m(D) 
m -m 
+n {1 = y ------ e ) 
x=0 ae 
If the supply of reserve items is very large, or large 


enough so that the chance of its exhaustion is negligible, 
then the random variable X is Poisson distributed with mean 
m(D). Therefore, 
CARR = ou pg es CLOG) amar ae 
iPiatse seamounts 
when the time permod. 15 80 ro) Db. 

For this scenario of item replacement, it is wsefiul ae 
look at CARF expressions for the various forms of intensity 
functions introduced earlier in this chapter. As an example 
of a linear intensity function, let a = 0, b = 0.04, D = en 


days and n = 25. Then 
ACE) = One4se 
and 


30 
m(30) = iP (0.04t)dt = 18 


Therefore, from (5.6) 


x De x 
neve! y oS -18 s 18 -18 
CARF = (---) ---- e + 25(1 - ---- e 
2S pe ame pg GL sed 
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= 12. Om 


Mins value can be approximated by (5.7) which results in 
fae Sete cae lemme edidmaere storm of Che intensity 


minmetion, let n = 20 and 
ety ea) 225 eee eee tt = - 1/225 t(t-30) 
Then 
30 
m(30) = ; (Ge 275 eee 2/15 edt = 20 
Therefore, the CARF = 91.12 4 which is obtained from (5.6). 
Mais value could also be approximated as 100 4% from (5.7). 


Examples of approximate and correct CARF values. are 


provided in Table 2. 


TAD LEZ. 


EXAMPLES SHOWING THE IMPACT ON CARF VALUES OF ASSUMING 
Pir ivGLe SPARES FOR A NONHOMOCENECGUS POISSON PROCESS 


CARF values(%) 


COvwoNns 
LININOUD 
WAU 


WdOnU, 
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Figure 5.5 suggests -that the approximation becomes 
better as the ratio m(D)/n becomes small. 


(D = 30 DAYS, N: NO OF IT 


Viz 
MEAN VALUES 





De Dns) The Impact _of Assuming 
nfinite Spares on CARF Values. 


So far, we have investigated ways of finding the CARF 
values for various scenarios according to several types of 
life distributions and for a nonhomogeneous Poisson process 
as an item's loss process. Based on this background, we 
will compare CARFs for each Jlife distribution and then 


Suggest some conclusions in the next chapter. 
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VI. COMPARISONS AND CONCLUSIONS 


oo - 


EweeecOMPARISONS OF CARFS FOR EACH LIFE DISTRIBUTION 


The purpose of comparisons between CARFs for each life 
PeewetDUtiOn is to find the type of life distribution which 
can be assumed for each combat scenario when we don't know 
mime exact life distribution of an item. Our decision rule 
mile be that, if we don't know the exact life distribution 
of an item, we should use the one yielding the maximum CARF 
in order to avoid the risk of underestimating the CARF 
value. All of the following figures were made by the values 


provided in the appendices. 


1. CARFs When Items Have the Same MTITL 


The first scenario we will examine is that when all 
items are on line with the same MTTL. 

Figure 6.1 shows CARFs as a function of MTTL values 
for the three distributions when all items have the same 
MTL. Note that these curves are simply plots of CARFs 
Seaemst MITLs and are not survival functions. The figure 
Suggests that the CARFs for the exponential are larger than 
those for the Weibull and gamma (shape parameter 2.0 and 
Be 0 ) when the MTTL is approximately greater than 20. 
Therefore, CARF values for the exponential could be used for 
Mets values ereater than 25 for this combat scenario. CARFs 
fore, tne Weibull distribution could be used for the MTTL 
values less than 25. 


2. CARFs Where Not All Items Have the Same MTTL 


mm rm 





When different subsets of the n items have different 


MTTL values, it is more difficult to reach general 


be 







(ALPHA=2.0, D=20 DAYS) (ALPHA=3.0, D=30 OAYS) 
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Figure 6.1 The Association Between CARFs and MTTLS 
When Items Have the Same MTTL. 


conclusions about which distribution it might be wusefuiagee 
assume. Our conclusions relates to the specific example 
examined earlier where there were three subsets of items, 
each With & ditterent MIlIL:. In Figure 6.2, the average Wor 
the three MTTL values was taken as an MTTL point. 

From Figure 6.2 for this scenario and example, we 
also realize that the CARFs for the exponential are the 
maximum values over the range of MTTL values greater than 
2% On the other hand, CARFs for the Weibull could be used 
for the MTTL values less than 25. 


ce CARFs Where There is a Change in the MITL 


a pe mm 


When all items have the same MTTL which changes once 
during the combat period, a partial comparison of distrisa. 
tions may be made using the numerical results developed 
earlier. In Figure 6.3, the average of the two MTTL values 


was taken as an MTTL point. 
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(ALERA=2.0, D=30 OAYS, P1=0.2, P2=0.5, F2=0.3) 


(ALSHS=3.9, 0=50 OAYS, F1=0.2, F2=0.8, F2=0.3) 
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Preure 6.2 The Association Between CARFs and MTTLs 
When Not All Items Have the Same MTTL. 
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There 1s a Change in the MTTL. 
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From Figure 6.3 for the particular example of tine 
scenario, we could see that the CARFs for the exponential 
are maximum values with MTITL values greater than 20 and 
CARFs for the Weibull are maximum with MTTL values less than 
Ebiss. 


B. CONCLUSIONS 


We have derived some mathematical expressions for CARFs 
in order to estimate equipment losses in the combat environ- 
ment by using the life distribution of a certain = ieee 
Furthermore, we discussed the way of finding the CARF value 
for the case where an item's loss follows a nonhomogeneous 
Poisson process. . 

We also attempted to suggest the kind of life distribu- 
tion that would be applied for each combat scenario by 
comparing these values when we don't know the exact life 
diStriputcion cisdmemtemn. 

CARF values that were drawn by this method could be uséd 
for procurement, SEOCKpalineg | and plans for shippiie 
requirements. However, the MTTL or the intensity function 
for an item must be estimated or obtained in advance to use 
this method. 

Consequently, using our conServative decision 7 rue 
seeking the larger CARF value, it appears from the computa- 
tions done here that the CARFs for the exponential instead 
of those for the Weibull or the gamma might be used if we 
don't know the exact life distribution of an item in which 
we are interested. However, for the smaller values of MTTL, 
CARF values for the Weibull appear to be better. 

In many cases, the population of n items which are 


considered for CARF valuation will be spread over variable 
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levels of combat intensity ,e.g., some items in the front 
and others in the rear while all are in the combat zone. It 
may be easier and more accurate to estimate MTTL values for 
each combat intensity level than to estimate an-MTTL for the 
population as a whole. Since the proportion of the n items 
at each intensity level may also be well estimated in many 
cases, the CARF expressions for the scenario where different 
subsets have the different MTTL values may be most useful. 
Mire choice of scenario and Jlife distribution will, ot 
course, depend primarily on the specific item being studied. 

The work on the preceding pages is an extension of past 
efforts on analytic modelling of CARF values by removing the 
need to assume constant loss rates or a homogeneous Poisson 
process. It is hoped that these extensions will be useful to 
those concerned with the accurate estimation of replacement 


feceors. 
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APPENDIX A 
CARFS WITH THE SAME MTTL (ALPHA = 2.0) 


(DO =s sOedays) 


MITitdays) CARF values(%) 

Exponential Weibull Gamma 

10.0 9520 Sie) ae 98.3 
20.0 ar a ase) 8. Ogi 
30R6 ba). 2 54.4 5 Sige 
40.0 52s Ce Ene 44.2 
50.0 45.1 24.6 33 
BOO Os Teg aes) 26.4 
fOs0 34.9 13.4 2 ee 
80.0 Sees, LOD lie 
90.0 aps eS 8.4 14.4 
LOO a0 pee, eva Lae 
LOO Za oye, 7 10.4 
jee e118) ZO: 4.8 an (0, 
5070 20.6 4.1 ia 
140.0 NES es: Shee) Ge 
1502 0 Loa t 6 ae Ge 
Lo@20 ae Co 5S 
EO O oer 2.4 4.9 
307 0 15.4 Ze 4 
EQORO 14.6 9 4.1 
20020 NES pee is Bie i 


_— = se eee os i i me ee eS se eee se ee eee ee eel ee 


44 


APPENDIX B 
SAR rome ine sAME MPTE (ALPHA = 3.0) 


(D = 30 days) 


MTTL(days) CARF values (%) 

Exponential Weibull Gamma 

LOG Ia 0 LOO. 99.4 
Z020 ims: 2 eo Seals) 
30.0 Soe 510 We Sie 
40.0 SS) aan ies ee) Sol 
5) ONO) 45.1 14.3 ILS ee, 
60R 0 BSS) S25 Toe 1 
BOrae 34.9 Dia) 14.0 
80.0 eile: 3g ff he => 
90.0 Zoro 256 80 
100.0 IL 2 eg Sle 
ARO © MB 1.4 D0) 
O30 Zo lage 4.1 
LS Or TAOS Ie. OR. SS, 
140.0 Eee: 0.7 Ze 
502.0 Ikeys I 0.6 Ze 
NGO70 ee 022 2,0 
Or. ieee 0.4 aay 
ESO 0 15.4 Dies 1.4 
SO) 14.6 O23 eee 
200.0 Lo OMe alee 
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APPENDIX C 
CARFS WHERE NOT ALL ITEMS HAVE THE SAME-MTTL 


(D = 30 dayemmpl= "90.2. vo2— 0. 52. eps Orem 


MTTL(days) | CARF values (%) 
MITLI °° Mia ye nie Exponential Weibull Gamma 


Te 20 S20 24.0 87.16 89.94 Sota 
52.0 Iter 18. 43.0 68.44 64.74 6 Gieze 
oa 0 S240 260 45.96 30s 364562 
23520 64.0 192.0 ET AS Feet 14.87 
256.0 2S 3Omeoc4. 0 ba 9 1 2a, 4.84 
(alpha = 3.0) 
IES Pe, 8.0 24.0 Gy Le D2 ce $3 ao 
SAD Ie a8 43.0 68.44 G27 72a 21.36 
64.0 5220 700.0 45.96 24.25 0.00 
Ze 6405 F920 Vine Ma BoC, 0.00 
2562.0 | 2 Se Onrmes 4. © ra I 49 0.00 
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APPENDIX D 
CARFS WHERE THERE IS A CHANGE IN THE MTTL 


wma O days, Dl = 15 days) 


MTTL(days) CARF values(%) 
Min 1 EE) 22 Exponential Weibull Gamma 


5,8 Lon Shee ei) 100.00 O96 5 
ORO ZOO 89.46 Se) oder she) 
2020 40.0 oes)! 53.84 62e56 
40.0 8020 43.02 py ees) a TANS, 
80.0 IoC 9) 24.52 Ae Fe. Oe 

160.0 Se 0 as 2 1220 je {oie 
(iboha = 3.0) 

0) 10.0 95.89 L0G 366 oo 31. 
LO 20 2020 89.46 vee ch9) Oia a7 
2.0.0 40.0 Seno) 43.07 Biche ge. 
40.0 OE 8; 43.02 6.80 18.24 
a0. 0 1h6 O20 24.52 0.88 ap oly, 

E60. 0 S726 219) See Oi: 0.56 


—~ ee es ew ew es em oe es oes ew ee i i i i i i 
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